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Abstract 

We consider operators with random potentials on graphs, such as the lattice 
version of the random Schrodinger operator. The main result is a general bound 
on the probabilities of simultaneous occurence of eigenvalues in specified distinct 
intervals, with the corresponding eigenfunctions being separately localized within 
prescribed regions. The bound generalizes the Wegner estimate on the density of 
states. The analysis proceeds through a new multiparameter spectral averaging 
principle. 
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1 Introduction 



In this note we present some general bounds for the joint distribution of eigenfunctions 
of operators with random potential, in the discrete setting. After finding that a natural 
multi-level extension of the Wegner bound on the density of states is not generally 
valid, we present a corrected version which is. It consists of a general bound on the 
probabilities of simultaneous occurrence of eigenvalues in distinct intervals, with the 
corresponding eigenfunctions being separately localized within prescribed regions, in 
a sense made precise below. The bound is derived through a suitable multi-parameter 
extension of the spectral averaging principle which is a familiar, and useful, element of 
the mathematical theory of Anderson localization. 

Wegner-type bounds are of relevance for the analysis of the extensions of the 
Schrodinger evolution to non-linear time evolutions and to interactive extensions of the 
one-particle model. These are not discussed here, but let us note that such systems con- 
tinue to attract attention, with interesting results presented in ref. IIAF88I IB W071 iFKSl 
as well as in a number of works which are currently in progress IIFKSllCSl . 

More explicitly, we consider random operators acting in the Hilbert space P{h), 
with A a finite set, of the form 

Hi,{uj)^T + V{u) (1.1) 

where T is a fixed hermitian operator and the randomness, represented by lo, enters 
only through a diagonal matrix V{uj) = dia.g{Vx{uj))x£A- Here w is a variable taking 
values in a probability space (ri,P). The joint distribution of {Vx{uj)}xeA induces a 
probability measure on the space of realizations MI'^L For convenience, and without 
loss of generality we identify O with this space, with given by the natural co- 

ordinates. Expectation values with respect to the probability measure P will be denoted 
by E. 

By default, it will subsequently be assumed here that the joint distribution of the 
potential variables satisfies the following regularity condition; 

Assumption R: For each site x £ A, the conditional probability distribution of Vx, 
conditioned on {Vy}y-tx, is absolutely continuous with respect to the Lebesgue 
measure, and its density (i.e. the corresponding Radon-Nikodym derivative) is 
uniformly bounded by a constant, p^o ■ 

Among the general results which are known for such random operators, and which 
have already played useful roles in the mathematical analysis of Anderson localization 
and of the corresponding spectral statistics, are: 

1. Spectral simplicity; With probability one Ha{uj) has only simple, i.e., non- 
degenerate, eigenvalues. (A proof which does not rely on the more involved 
Minami estimate is spelled in the appendix.) 

2. The Wegner bound llWeSlI : the mean density of states is bounded by poo- 
Equivalently, for any energy interval / 

P{a(i/A)n/^0} < E[TtPi{Ha)] < Poo|/||A|, (1.2) 
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where a is the spectrum of the operator, Pj is the corresponding spectral projec- 
tion, and I • I denotes a set's length, or 'volume', as appropriate. 

3. The Minami bound ||Mi96| : the probability of there being multiple eigenvalues 
in a small energy range satisfies 

P ({card{cr(HA) n 1} > 2}) < ^[Tv Pi{H) {Tv Pi{H) - I)] 

< ypLl/PlAp. (1.3) 
(The statement had a one dimensional precursor in IMoSlll .) 
These bounds were recently extended IIGV071 IBHSOTI to: 

P ({card{a(HA) n 1} > n}) < J 1/1" lAI" ■ (1-4) 

Furthermore, in a work which was posted at the time of completion of this article, 
the Minami bound was given a new and more transparent derivation and some further 
extensions IICGK08I . 

At first glance, one could ask whether (11.4b is a special case of a more general valid 
bound on the n-point density functions, of the form: 

k 

P({a(i/A)n/, ^0 foralli = l,...,A:}) < C„p^ n|/,||A|, (1.5) 

where {Ij} could be arbitrary intervals. 

A bound like il.5\ could be of use, e.g., in estimating the probability that for an 
a-priori specified energy E there is a of multi-state resonance, in the sense that the 
quantity \E — ''^j-^jl is small for some integer sequences {rrij}. Questions of 

this kind are of relevance for the non-linear extension of the Schrodinger evolution 
which is studied in IFKSII . 

As it turns out, the bound suggested in dl.Sb does not hold at the generality of the 
two preceding statements. As the simple calculation which is demonstrated next shows, 
it is already not valid for the 2x2 example. Nevertheless, a somewhat similar bound 
does hold for disjoint energy intervals under the restriction that the eigenfunctions' 
moduli are of sufficiently different profiles. The precise statement, which is our main 
result, is presented in section[3] 



2 A counterexample 

While at first glance ( 11.5b may appear sensible, and even supported by the observation 
that for random matrices the level interaction is repulsive, it is easily seen to be false. 
A counterexample to ( 11.51 ) is found already in the context of 2 x 2 matrices. 

In the two dimensional space, a self adjoint operator with random potential is given 
by a 2 X 2 self adjoint matrix of the form: 

J. (2.0 
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with some a, & G R and c G C. The two eigenvalues of H{uj) are 

Ei/2{uj) = ^(uji+LU2 + a + b± V(wi-t^2+a-6)2+4|c|2) , (2.2) 
where one may note that the spectral gap satisfies: \Ei{uj) — E2{u!)\ > 2\c\ for all 

The determinant of the change of variables (0^1,^2) (£'i,i?2)is given by 



dct 



dEj{uj 
dtjJi, 



\oji~oJ2 + a-b\ _ ^{Ei- E2Y -^c\ 
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^{L0i-W2+a-bY+A\c\^ \Ei-E2\ 

(2.3) 



with i,k G {1, 2}. Hence, for lui and W2 a pair of iid variables with a common density 
Q, the probability density for the pair of eigenvalues £"2) (with Ei ^ E2) is: 



p{Ei,E2) 



det 



dujk 



q{0Jj{EuE2)) (2.4) 



0/ \E^-E2\<2\c\. 



with (i?! , E2) determined by the relation i 

The above density has the singularity of {\Ei — E2\ — 2|c|)^^/^ at the edge of 
the spectral gap, where \Ei — i?2| = 2|c|. In effect, one could see here that the level 
repulsion is associated with a 'pile - up' of the probability density at the edge of the 
gap which it creates. 

Clearly, similar singularities in the two point function would be found in the more 
general n x n situation whenever the systems is decomposable with an isolated two-site 
subsystem. However, it seems to be an interesting question whether for generic n x n 
matrices with random potential the singularity is rounded off due to the larger number 
of random variables. 

One may also note that while the above calculation contradicts dl.St . it implies that 
at least in the 2x2 case the two point function satisfies a modified bound, which is be 
obtained by replacing \Ij\, on the right side of ( 11. 5t . with max{|/j |^/^, \ It will be 
of interest to clarify how far can such a bound be extended. Suitable generalizations 
could provide useful information on the probabilities of multi-level resonances which 
are mentioned above. 



3 Bounds on the joint distribution of the eigenvalues 
3.1 Statement of the main result 

The 'positive results' which are presented here amount to bounds on the probabilities of 
simultaneous occurrences of eigenvalues, in prescribed intervals, which are associated 
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with eigenfunctions of sufficiently distinct profiles. Following is the definition of that 
concept. 

Definition 3.1. Normalized functions ■01 , . . . , ■(/'„ G with = 1, are said to 

have a-distinct profiles, for some a > 0, within sets . . . , B„ C A if and only if 



>a". (3.1) 



It may be noted that by the linearity of the determinant and the triangle inequality: 



(3.2) 



where Id stands for the indicator function of the set D. Hence, a sufficient condition 
for ( 13. Il l is that the row (or column) vectors in the (substochastic) matrix of occupation 
probabilities Ib^ k=i ^P^*^ ^ parallelepiped of volume at least a". 

We shall now consider the events: 

8a{Il- ■ ■ ■ , In', Bi, . . . , Bn) (3.3) 

II/^{lo) has eigenvalues E'l e /i, . . . , En € In whose eigen- 
functions have a-distinct profiles, within sets Bi, . . . , Bn, 

with /i , . . . , /„ C R a collection of Borel sets, and Bi, . . . , Bn C A. 

Proven below is the following statement. 

Tlieorem 3.1. For operators with random potential, as in (II. lb . whose probability 
distribution satisfies the regularity assumption R, the probabilities of the events defined 
in ( |3.3l l satisfy: 



'(£„(/!...., /„;i?i,...,B„)) < — (3.4) 



Concerning the uses of this result, it may be noted that throughout the localization 
regime, where the eigenfunctions are each localized in some region of -roughly- the lo- 
calization length, condition ( 13.1b would be satisfied by eigenfunctions of separate sup- 
ports. Theorem 13. II can therefore by used to bound the probabilities of eigenfunctions 
in prescribed intervals whose eigenfunctions do not overlap in space. Regrettably, the 
results presented here do not address the corresponding question for eigenvalues with 
overlapping eigenfunctions. One may wonder whether even in that case the eigenfunc- 
tions' profiles should typically be distinguishable, in the sense of ( 13.1b . Such a result 
could extend the applicability of Theorem l3.1l 



3.2 Multiparameter spectral averaging 

To prove Theorem |3.1| we first establish the following estimate. 
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Lemma 3.2 (Generalized spectral averaging). For operators with random potential, as 
in d 1 . 1 b , whose probability distribution satisfies the assumption R, and any collection 
of intervals /i , . . . , /„ and sites xi^ . . . ,x„ G A; 



E 



det((5,,,P/^(ffA)<5.,)) 



< nl 



n i^ii 



(3.5) 



The term by which we refer to this statement is motivated by the observation that 
the case n = 1 yields the known spectral averaging principle: 



(3.6) 



from which the Wegner estimate ( 11.2b readily follows. The proof of the more general 
statement is based on an elementary change of variable calculation, combined with in- 
put from algebraic-geometry. The latter is needed for a bound (which is rather natural) 
on a relevant multiplicity factor 

Proof of Lemma [O] We shall first derive (13.5b under the additional restriction to the 
event: 



Jih 



The spectrum of H\{llj) includes exactly one 
eigenvalue in each of the intervals /i, .., /„ 



For u! e ,...,/„) the determinant in the left side of ( 13. 5b reduces to 

D(E;Y) det„x«(|V'j(a;fc)n , S := {xi, . . . ,a;„} , 



(3.7) 



(3.8) 



where E ;= (i^i, . . . .En) is the set of eigenvalues which occur in the indicated inter- 
vals, and ijjj are the normalized eigenvectors of H\{lu) corresponding to the (uniquely 
defined) eigenvalues Ej E Ij. Thus, our first goal is to establish the bound 



E[\D{E;J:)\ 1j,,„.„,J < 



(3.9) 



where Ij = Ij-(w) denotes the indicator function of the event J" = . . . , /„). 

The expectation value in ( 13. 9b can be calculated as the average of the conditional 
expectation of the same quantity conditioned on Va\s {Vx}x^'s, i.e. 



E[\D{E;mij] ^ f ( Ij \D{E-Y.)\ ii{dV^\V!,\^) 

Thus, under the assumption R on the joint distribution of {T4}, we have: 
E[\D{E-Y)\lj] < sup / \D{E;j:)\dV^. 

where S is the following subset of the section of J at the specified Va\s 
S:={Vj:\v= (Fs, 14\s) e J and D{E; S) ^ } . 



(3.10) 
(3.11) 

(3.12) 
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The integral on the right side of ( 13.1 11 1 may be conveniently expressed through the 
change of variables 



■.^{V^„...,V^J^E:={Ei,...,En), (3.13) 

which is to be understood as performed at fixed Va\s. Standard perturbation theory 
IIKa66l implies that the set J' C M''^' is covered by open sets within each of which 
Ej, are defined as single-valued analytic functions of Vs, with derivatives given by the 
Feynman-Hellmann formula: 

f) E 

' - m^k)\' . (3.14) 



Hence, the Jacobian for the coordinate change is given by 

The section S is covered by open sets on which D(^E; S) ^ 0, for which the transfor- 
mation ( 13.131 ) is locally bijective. Globally, the mapping is not 1 — 1, and the correct 
change of variables formula is; 

[ \D{E;J:)\ dVj: = [ N{E:Y.) dE . (3.16) 

with the multiplicity factor: 

N{E; S) := card {Vs | E are eigenvalues of Ha{Vy.,Va\s) and D{E; S) 7^ O} . 

(3.17) 

The factor N{E; S) counts the number of simultaneous solutions, for Vs, of the 
set of equations (at fixed Va\s)- 

PeAVx^,-.-,V^J = 0, j = l,...,7i, (3.18) 

where Pe{Vxi , ■ ■ ■ , K-„ ) is the characteristic polynomial: 

P£;(14i,...,14J = P(14,,...,14„;S) dct{HA-E) (3.19) 

The number of solutions of a system of algebraic equations is a classical problem of 
algebraic geometry (for which it is the size of a zero dimensional algebraic variety, 
defined by ( 13.181 )). A rather general answer is provided by so-called Bezout's theory. 
However, a simplifying observation is that in the case of interest for us the polyno- 
mials Pe{Vxi , . . . , Vx„) are linear in each of the Vxj variables. To form a guess as 
to what may be the number of solutions for such systems, one may observe that the 
answer is trivial if the non-random term in H/^{uj) = T + V{uj) is a diagonal matrix 
diag{Txi, ...,Tx„}. In this case, equations ( 13.18b are simultaneously satisfied if and 
only if {Tx + }j=i,....n coincide with a permutation of {£'j}j=i....,n, and thus 
N{E;Y,) <n\. 
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As it turns out, by a theorem due to D. Bernsteird also in the more general case 
which is of interest to us the number of solutions of the system (13.18b satisfies; 

N{E;J:) < nl . (3.20) 

The applicable theorem is Proposition ! A. 1 1 which is presented in AppendixlAl To apply 
it, we need to check that for the counted solutions det ^ ^gy'^' j 7^ 0. For that we note: 

(3.21) 

Since the last product is non-zero for ?i\\V £ J and j = 1 , . . . , n, condition ( IA.3I ) is 
satisfied on S. 

The above considerations prove (I3.9l l, which differs from ( I3.5l l mainly in the pres- 
ence of the additional constraint that each interval Ij includes exactly one eigenvalue 
of Ha. We shall now show that ( 13.5b follows. For that, consider a partition of Uj/j into 
a finite collection, C^, of disjoint sub intervals whose length does not exceed e. One 

I ■ (e) 

may represent each of the intervals Ij as a disjoint union Ij = u;;,^!/™ of elements of 
such a partition, i.e. with each in (if the sets Ij are not pairwise disjoint some 
elements of will be called upon more than once.) By the linearity of the determinant: 

det ((4. , Pi, {Ha) <5., )) = ^ ... ^ det ((<5,, , P^^ (Ha) 4. )) , (3.22) 

mi — 1 m^i — 1 

The sum can be restricted to the case that , . . . , In" are disjoint, since otherwise 
the determinant vanishes. 

Using the fact that the determinant on the left is bounded by one, we can estimate: 

In at least one of the elements of 
Ha has two or more eigenvalues 



E[|det((4.,P7,(i?A)4.))|] < 



(3.23) 



where the summation range is as in ( 13.221 ). We shall now take the Umit e ^ 0. Since 
the eigenvalues of Ha are almost surely simple (Lemma lB.lb . the dominated conver- 
gence theorem implies that in that limit the first term vanishes . Applying (13.9b to the 
remaining terms one gets: 

E (I det ((4. , Pi, (Ha) 4. )) |) < 1™ ^ ... E IT 1^^' I 

mi — 1 m^ — 1 j 

^n\p2,ll\I,\ (3.24) 



^We thank J. Kollar for help with the reference. 
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which yields ( 13.5b . 



□ 



3.3 Proof of main result 

Proof of Theorem \3.1\ Using the limiting argument employed at the end of the proof 
of Lemma 13.21 one shows that it is sufficient to bound the probability of n J, 
where J was defined in ( 13.7b . Using the assumption on the normalized eigenfunctions 

ipi, . . . ,ipn corresponding to Ei e /i, . . . , £"„ G In,we then estimate: 

P ...,/„; Bi, ... , S„) n J{h. /„)) 

E E lE[|dct(|^,(xfc)|2)| l^(,,.....z„)] 

71 

<a-n\Ql,\{\I,\\B,\ (3.25) 

i=i 

where the last inequality is due to (13.9b . □ 



Appendix 

A Counting solutions of a system of polynomial equa- 
tions 

In the proof of Lemma 13.21 for the bound (13.201) we invoked the following general 
result. 

Proposition A.l (Special case of a theorem by D. Bernstein). Let Pj, j = 1, ...,n 

be polynomials in n variables, cr = (ai, • • • ,cr„) G C", which are linear in each 
variable, i.e., are of the fonn 

P,{a)= J2 , (A.l) 

fc6{04}" 

with Cj{k) e C which are non-zero only if km ~ 0, 1. Then the number of isolated 
solutions of the system 

Pj{(t)^0, forallje{l,...,n}, (A.2) 
is at most nl. Moreover, each solution cr at which 

dp,ia)Y 



det ^0, (A.3) 
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is isolated. 



The first part is a special case of IICLS91I Thm. 5.4]. In case the solution is not 
isolated, there exists locally a differentiable curve s i-^ c(s) such that for all j = 
l,...,n 

rfPj-K,s)) _^ dP,{ais))dakis) ^ ^^^^ 

ds ^ d(Ti. ds 

fc=i " 



This contradicts the assumption ( IA.3l l. which implies that the matrix of partial deriva- 
tives has no zero eigenvalue. 



B Simplicity of the spectrum 

In our discussion it was convenient to know that the spectrum of an operator with 
random potential is almost surely non-degenerate. While this assertion is among the 
consequence of the Minami bound, for completeness we present here also a direct and 
elementary proof. 

Lemma B.l. Let H\(uj) be an operator in ^(A)^, for some finite region |A|, with 
a random potential such that for each x € A the conditional distribution of Vx{uj), 
conditioned on {K^(w)}.yg|^}c, is almost surely continuous. Then for almost all uj the 
spectrum of H\{uj) has only simple eigenvalues. 

Proof of Lemma \BJ] Let ipi, . . . , tp\j^^ be an orthonormal basis of eigenfunctions of 
Ha with corresponding eigenvalues denoted hy Ei, . . . , E^\\. Consider the self adjoint 
operator 

A/a := (i/A ® 1 - 1 ® Ha)^ 

on the subspace 7i ~ of antisymmetric functions within the product space €^ (A)®^^ (A). 
It is straightforward to check that the orthonormal basis given by 

*7fc ^ (V-j ® V^fc -i'k® i'j) , .] <k, (B.l) 

constitutes an eigenbasis with A/a ^'Jj. = [Ej — Ek)'^ '^Jf^. The simplicity of the 
spectrum of Ha is therefore equivalent to A/a being almost surely invertible on , 
i.e., 

dct A/a > . (B.2) 
For a proof of this assertion, we consider the ' 2 ' ^ matrix-elements given by 

{S:,y,,MAS-y), 5-y 1= {5^ (^i 6 y - 6 y (^i S (B.3) 

associated with the antisymmetrized position basis of 7i~, and study their dependence 
on a single random variable, say Vx- Only |A| — 1 rows (and columns) of the matrix 
depend on 14. In these rows, the diagonal matrix elements, with x ^ x' and y ~ y', 
depend on Vx quadratically, while the diagonal elements are linear in Vx- 
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Hence, Vx i— *■ detMA is a polynomial of degree at most |A|. Thus, we have the 
following dichotomy: the characteristic polynomial has either at most | A| isolated zeros 
or is independent of 14. It the first case, the conditional probability that dot AIj\ = 
0, conditioned on {Vy}y^x, vanishes, since the distribution of Vx is assumed to be 
continuous with respect to Lebesgue measure. In the second case, one may reduce the 
site X from A by taking the limit Vx oo. In this limit Ha ^ H\\^x} © oo. Since 
dot M\ does not diverge in this limit, one may conclude that dct Ma\{o} = 0, and the 
argument may be repeated for the smaller set. In case |A| = 1 the condition iB.2i is 
trivially satisfied. This completes the proof. □ 
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